Coordination of passive systems 
under quantized measurements * 

Claudio De Persis ^ Bayu Jayawardhana ^ 
April 10, 2012 

Abstract 

In this paper we investigate a passivity approach to collective co- 
ordination and synchronization problems in the presence of quantized 
measurements and show that coordination tasks can be achieved in a 
practical sense for a large class of passive systems. 

1 Introduction 

In the very active area of consensus, synchronization and coordinated control 
there has been an increasing interest in the use of quantized measurements 
and control ([28, 31, 24, 8, 29, 11] and references therein). As a matter of 
fact, since these problems investigate systems or agents which are distributed 
over a network, it is very likely that the agents must exchange information 
over a digital communication channel and quantization is one of the basic 
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limitations induced by finite bandwidth channels. To cope with this lim- 
itation, measurements are processed by quantizers, i.e. discontinuous maps 
taking values in a discrete or finite set. Another reason to consider quantized 
measurements stems from the use of coarse sensors. 

The use of quantized measurements induces a partition of the space of mea- 
surements: whenever the measurement function crosses the boundary be- 
tween two adjacent sets of the partition, a new value is broadcast through the 
channel. As a consequence, when the networked system under consideration 
evolves in continuous time, as it is often the case with e.g. problems of coor- 
dinated motion, the use of quantized measurements results in a completely 
asynchronous exchange of information among the agents of the network. De- 
spite the asynchronous information exchange and the use of a discrete set of 
information values, meaningful examples of synchronization or coordination 
can be obtained ([20, 13, 30, 18]). 

There are other approaches to reach synchronization (consensus) with an 
asynchronous exchange of information, such as gossiping ([7, 9]), where at 
each time step two randomly chosen agents exchange information. The lat- 
ter approach is conceptually very different from the one considered in this 
paper, where the agents broadcast information when a local event occurs 
(the measurement crosses the partition boundary). Moreover, while gossip- 
ing algorithms are mainly devised for discrete-time systems, here we focus 
on continuous-time systems. 

In view of the several contributions to quantized coordination problems 
available for discrete-time systems ([28, 31, 24, 8, 29, 11]), one may won- 
der whether it would be more convenient simply to derive the sampled-data 
model of the system and then apply the discrete-time results. Due to the 
distributed nature of the system, a sampled-data approach to the design of 
coordinated motion algorithms presents a few drawbacks: it might require 
synchronous sampling at all the nodes of the network and consequent accu- 
rate synchronization of all the node clocks; it might also require fast sampling 
rates, which may not be feasible in a networked system with a large number 
of nodes and connections. Finally, the sampled-data model may not fully 
preserve some of the features of the original model. For these reasons, we 
focus here on continuous-time coordination problems under quantized mea- 
surements. 

A few works on this class of problems have recently appeared. The work [18] 
deals with consensus algorithms using binary control algorithms. In [20] the 
attention is turned to quantized measurements and the consensus problem 
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under quantized relative measurements is tackled. The same problem, but 
considering quantized absolute measurements, is studied in [13]. The paper 
also introduces hysteretic quantizers to prevent the occurrence of chattering 
due to the presence of sliding modes. More recently, the work [30] has stud- 
ied the quantized consensus algorithm for double integrators. A remarkable 
advancement in the study of consensus algorithms over time-varying com- 
munication graphs and using quantized measurements has been provided by 
[23]. 

Despite the unquestionable interest of the results in papers such as ([20, 
13, 30, 18, 23]), they present an important limitation: they focus on agents 
with simple dynamics such as single ([20, 13, 18, 23]) or double integrators 
([30]). The goal of this paper is to investigate the potentials of an approach 
to coordinated motion and synchronization which takes into account simul- 
taneously complex dynamics for the agents of the network and quantized 
measurements. 

In coordinated motion, variables of interest are the position and the velocity 
of each subsystem, and the problem is to devise control laws which guarantee 
prescribed inter-agent positions and velocity tracking. In this paper we focus 
on the approach to coordinated motion proposed in [1]. In that paper, the 
author has shown how a number of coordination tasks could be achieved for a 
class of passive nonlinear systems and has been using this approach for related 
problems in subsequent work ([5, 6]). Others have been exploiting passivity 
([14, 26, 42, 32] to name a few) in connection with coordination problems. 
Our interest for the approach in [1] stems from the fact that it allows to deal 
with complex coordination tasks, including consensus with velocity tracking, 
in way that naturally lends itself to deal with the presence of quantized mea- 
surements. In the approach of [1], a continuous feedback law is designed to 
achieve the desired coordination task under appropriate conditions. Thus 
the presence of quantized measurements can be taken into account in this 
setting by introducing in the feedback law static discontinuous maps (the 
previously recalled quantizers). Although in the case of quantized measure- 
ments the conditions in [1] are not fulfilled due to the discontinuous nature 
of the quantizers, one can argue that an approximate or "practical" ([13]) 
coordination task is achievable under suitably modified conditions. This is 
the idea which is pursued in this paper. In the case of a control system with 
a single communication channel this was studied in [12]. Another reason to 
consider the approach of [1] is that it provides a systematic way to deal with 
a large variety of cooperative control problems, as it has been authoritatively 
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proven in the recent book [3]. 

A second aim of this paper is to study practical state synchronization under 
quantized output feedback. In these problems, one investigates conditions 
under which the state variables of all the subsystems asymptotically con- 
verge to each other, with no additional requirement on the velocity tracking. 
Passivity ([14, 39, 38]), or the weaker notion of semi-passivity ([36, 35, 40]), 
has also played an important role in synchronization problems. Here we 
mainly focus on the models considered in [14, 38]. 

The main contribution of this paper is to show that some of the results of [1] 
and [38] hold in a practical sense in the presence of quantized measurements. 
Because the latter introduces discontinuities in the system, a rigorous analy- 
sis is carried out relying on notions and tools from nonsmooth control theory 
and differential inclusions. As far as the coordination problem is concerned, 
although the passivity approach of [1] allows to consider a large variety of 
coordination control problems, in this paper we mainly focus on agreement 
problems in which agents aim at converging to the same position. 
A few other papers have appeared which deal with coordination problems for 
passive systems in the presence of quantization. The work [26] deals with a 
position coordination problem for Lagrangian systems when delays and lim- 
ited data rates are affecting the system. The paper [22] deals with master- 
slave synchronization of passifiable Lurie systems when the master and the 
slave communicate over a limited data rate channel. The main difference of 
our paper compared with [26, 22] is that in the former each system in the 
network transmits quantized information in a completely asynchronous fash- 
ion and no common sampling time is required. Prom a mathematical point of 
view, this means that our approach yields a discontinuous closed-loop system 
as opposed to a sampled-data one. Moreover, the classes of systems and the 
coordination problems considered here appear to be different from those in 
[26, 22]. 

The organization of the paper is as follows. The passivity approach to co- 
ordination problems is recalled in Section 2. In Section 3 the coordination 
control problem in the presence of uniform quantizers is formulated and the 
main results are presented along with some examples. The synchronization 
problem for passive systems under quantized output feedback is studied in 
Section 4. In Section 5 a few guidelines for future research are discussed. 
In the Appendix some technical tools are reviewed for the sake of readers' 
convenience. 
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2 Preliminaries 



Consider N systems connected over an undirected graph G — {V, E), where 
y is a set of N nodes and E ^ V x V is a. set of M edges connecting the 
nodes. The standing assumption throughout the paper is that the graph G is 
connected. Each system i, with i = 1, 2, . . . , A^, is associated to the node i of 
the graph and the edges connect the nodes or systems which communicate. 
Each system i is described by 

where the state G M"% the input Ui G the output wi G W, the 
exogenous signal Vi G W and the maps fi,gi,hi are assumed to be locally 
Lipschitz satisfying /j(0) = 0, gi{0) full column-rank, /ij(0) = 0. For the 
system Ej, we assume the following: 

Assumption 1 There exists a continuously diff'erentiable storage function 
Si : M"'' M+ which is positive definite and radially unbounded such that 

^SM^){M^^) + gi{^i)Ui) < -WiiCi) + hiiCifu^, (2) 

where Wi is a continous positive function which is zero at the origin. 

Such a system Sj is called a strictly-passive system (with Vi — 0). If Wi is a 
non-negativefunction, then Ej is called a passive system. 
Label one end of each edge in £^ by a positive sign and the other one by a 
negative sign. Now, consider the k-th edge in E, with k G {1,2, .. . ,M}, 
and let i,j be the two nodes connected by the edge. For the coordination 
problem, which is detailed in Subsection 2.1, the relative measurements of 
the integral form ■u;j(r)dr and Wj(r)dr are used. On the other hand, for 
the synchronization problem, which is briefly reviewed in Subsection 2.2, we 
need the relative measurements of the signals Wi and wj. Thus, depending 
upon specific problems, let Zk describe the difference between the signals Wi 
and Wj (or the difference between the signals Xi{t) :— Wi{T)dT + Xi{0) and 
Xj{t) := J*Wj{T)dT + Xj{0) with constant vectors Xi{0),Xj{0) G W) and be 
defined as follows: 

{Wi — Wj (or Xi — Xj ) if i is the positive end of the edge k 
Wj — Wi (or Xj — Xi) if i is the negative end of the edge k . 
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Recall also that the incidence matrix D associated with the graph G is the 
N X M matrix such that 



+1 if node i is the positive end of edge k 
dik = { ~1 if node i is the negative end of edge k 
otherwise. 

By the definition of D, the variables z can be concisely represented as 

z = (D^ ® Ip)w (^or z = {D'^ ® Ip)x ^ (3) 

where w = [wf . . . wjf]'^ and x = [xf . . . xj^Y ? respectively, and the symbol ® 
denotes the Kronccker product of matrices (see Appendix A for a definition). 
In this paper we are interested in control laws which use quantized measure- 
ments. For each A; = 1, 2, . . . , M, instead of Zk, the vector 

q(^fe) {q.{zki) ■ ■ ■ q.izkp))'^ 

is available, where q is the quantizer map which is defined as follows. Given 
a positive real number A, we let q : R — >■ ZA be the function 



q(r) = A 



r 1 
A + 2 



(4) 



with ^ the precision of the quantizer. As A — > 0, q(r) — )■ r. Observe that 
each entry of z^ is quantized independently of the others and the quantized 
information is then used in the control law. 



Remark 1 The results of the paper continue to hold if each quantizer has 
its own resolution (that is, the information Zkj is quantized by a quantizer 
with resolution Akj). However, to reduce the notational burden, we only deal 
with the case in which the quantizers have all the same resolution A. 

In the following subsections, we review the results on passivity approach 
to the coordination problems of [1] and to the synchronization problems of 
[38] without the quantized measurements. 
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2.1 Passivity approach to the coordination problem 

In the coordination problems of [1], the signal Wi of each system Ej corre- 
sponds to the velocity of the system, and thus, Xi, i — 1, . . . , N, represents the 
positions which must be coordinated (recall that Xi(t) := Wj(r)dr + ,Xj(0)). 
The coordination problem under consideration requires all the systems of the 
formation to move with a prescribed velocity v, i.e., vi = V2 = ■ ■ ■ = vn = v. 
Define 

y.^Xi-V (5) 

the velocity tracking error. It can be checked from (1) and the definition 
of Xi that Hi = /i($t). The standing assumption is that, possibly after a 
preliminary feedback which uses information available locally, each system 
Ej is strictly passive, i.e., (2) holds with Wi positive definite. In other words, 
it is strictly passive from the control input Ui to the velocity error yj. 
For the sake of conciseness, the equations (1), (5) are rewritten as 

/ /ii(ei) \ / v\ 



X 



+ 



~~m' 

( mi) 



\ 







(6) 



+ 







u 



9N{iN) I 



/(?) 



9{0 



[X^ 



X 



N 



TlT 
Ni ; 



U 



[ui ...ujf] , Iat is the A^- 



where x 

dimensional vector whose entries are all equal to 1 and denotes a vector of 
appropriate dimension of all zeros. 

The formation control problem consists of designing each control law Ui, 
with i = 1, 2, . . . , A^, in such a way that it uses only the information available 
to the agent i and guarantees the following two specifications: 

(i) lim(_>oo \xi(t) — v(t)\ = for each i = 1,2, . . . , N, with v{t) a bounded 

and piece- wise continuous reference velocity for the formation; 

(ii) Zk{t) Ak as t — oo for each k = 1,2,..., M, where Ak C MP are the 

prescribed sets of convergence ^ and z = {D^ (E) Ip)x as defined in (3). 



^We refer the interested reader to [1] for examples of sets Ak related to some coordi- 
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In [1], where measurements without quantization are considered, the case 
Ak = {0} is referred to as the agreement problem. 

Let Pfc : — )■ M, for k = 1,2, . . . , N , be nonnegative continuously dif- 
ferentiable (the latter assumption will be removed in the next section) and 
radially unbounded functions whose minimum is achieved at the points in 
Ak- To be more precise, the functions Pk are assumed to satisfy 

Pk{zk) = and V Pk{zk) = if and only if Zk G Ak- (7) 

Define 

WPk{zk) = ipk{zk) - (8) 
The feedback laws proposed in [1] to solve the problem formulated above are: 

M 

-^dikipk{zk) , 1^1,2,..., N. (9) 
fe=i 

Observe that, as required, each control law Ui uses only information which is 
available to the agent i. Indeed, dik 7^ if and only if the edge k connects i 
to one of its neighbors. In compact form, (9) can be rewritten as 

u ^ -{D ^ Ip)7/j{z) , (10) 

where iIj{z) = [ipilzi)'^ . . .i/jMi^M)'^]'^ and z is as in (3). Before ending the 
section, we recall that the system below with input x and output —u, namely 
(see Figure 2 in [1] for a pictorial representation of the system) 

z = {D'^(g)Ip)x . s 

-u = {D®Ip)i>{z) ^ ' 

is passive from x to —u with storage function X^^^ Pk{zk)- We remark that 
the function Pk{zk) is chosen in such a way that the region where the variable 
Zk must converge for the system to achieve the prescribed coordination task 
coincides with the set of the global minima of Pk{zk)- Hence, the coordination 
task guides the design of Pk{zk) which in turn allows to determine the control 
functions (9) via (8). The functions Pk{zk) in the case of agreement problems 
via quantized control laws will be designed in Section 3. 

nation problems. The sets Ak which are of interest in this paper will be introduced in 
(16). 
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2.2 Passivity approach to the synchronization problem 



In the synchronization problem of [38, Theorem 4], each system Ej in (1) 
(with Vi — 0) is assumed to be hnear, identical and passive. For such setting, 
each (passive) system Ej is of the form 

Wi = C6 i = l,2,...,7V ^ ^ 

where G W\ Ui,Wi G and the passivity of Ej implies that the following 
assumption holds: 

Assumption 2 There exists an {n x n) matrix P — > such that 

A^P + PA < 0, B'^P = C . 

The synchronization problems can then be stated as designing each con- 
trol law Mj, i = 1, 2, . . . , iV, using only the information available to te agent 
i such that, for every — ^ A where .^o is the trajectory of the au- 
tonomous system — A^q which is initialized by the average of the initial 
states, i.e., ^o(O) = I^i6(0), and ^ C is the prescribed set of conver- 
gence. In the case without the quantized measurements, which is treated 
in [38], A = {0}. The coordination problem that is reviewed in Subsection 
2.1, is related to the case when = [1]. For another viewpoint, we can 
consider that (12) corresponds to the case in the Subsection 2.1, where the 
mapping u y is an identity operator, v — and one takes into account 
dynamics on the subsystem x which are more complex than those of a single 
integrator. 

In addition to output synchronization, it is well-known that the states 
of interconnected passive systems synchronize under observability assump- 
tion ([14]). The largest invariant set of the interconnected systems, when 
the measurements are not quantized and (C, A) is observable, is the set 

G M"^ : ^1 = . . . = ^n}- In the case of quantized measurements, the 
invariant set is larger. Our main result in Section 4 provides an estimate 
of the invariant set of the interconnected systems with quantized measure- 
ments. To this purpose, we rely on a result of exponential synchronization 
under static output feedback control laws and time- varying graphs which has 
been investigated in [38]. In the following statement, we recall Theorem 4 of 
[38] specialized to the case of time-invariant undirected graphs: 
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Theorem 1 Let Assumption 2 hold and suppose that the pair (C, ^4) is ob- 
servable. Let the communication graph be undirected and connected, and 
denote z = {D^ ® Ip)w as in (3) with w = [wf . . . wjf]'^. Then the solutions 
of 



M 



i,^Aii-BY,(hkZk, i = l,2,...,7V (13) 

k=l 

satisfy 



lim 

t->-+oo 



(14) 



where ^ — [^1^2 ' ' ' ^JrV ^'^^ convergence is exponential. More precisely, 
the solutions converge exponentially to the solution of = A^q initialized to 
the average of the initial conditions of the systems (13), i.e. ^o(O) — 1% <^ 

inm/N. 

Let I = C - l^l^M^^ = (n ® with n = /at - be the disagree- 

ment vector. From (13), $,{t) obeys the equation 

i = [lN®A-(lN®B)(DD'^®Ip)(lN^C)]i 

A 

and the convergence result (14) can be restated as lim^^+oo ||C(^)|| = 0. The 
proof of the result rests on showing that the Lyapunov function 

along the solutions of (15) satisfies the inequality 

T>(0<-A2||(n®7,)|f, 

where A2 is the algebraic connectivity of the graph, i.e. the smallest non-zero 
eigenvalue of the Laplacian L — DD^. Then the thesis descends from the 
observability assumption and Theorem 1.5.2 in [37]. 



3 Quantized coordination control 

3.1 A practical agreement problem 

Despite the generality allowed by the passivity approach of [1] , in this paper 
we focus on an agreement problem. By an agreement problem it is meant 
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a special case of coordination in which all the variables Xi connected by a 
path converge to each other. In the problem formulation in Section 2, this 
amounts to have Ak = {0} for all /c = 1, 2, ... , M. When using statically^ 
quantized measurements, however, it is a well established fact ([28, 20, 13]) 
that a coordination algorithm leads to a practical agreement result, meaning 
that each variable zj^ converges to a compact set containingthe origin, rather 
than to the origin itself. Motivated by this observation, we set in this paper a 
weaker convergence goal, namely for each A; = 1, 2, ... , M, we ask the target 
set Ak to be of the form: 

Ak= X [-a, a] (16) 

where a is a positive constant and the symbol x denotes the Cartesian prod- 
uct. Then the design procedure of Section 2 prescribes to choose a non- 
negative potential function Pk{zk) which is radially unbounded on its domain 
of definition and such that (7) holds. If such a function exists then the con- 
trol law is chosen via (8). To take into account the presence of quantized 
measurements, the nonlinearities on the right-hand side of (8) should take 
the form 

^k{zk) = Xik(q(^fc)) , (17) 

with Xk to be defined later. 

The presence of quantized measurements, i.e. of q(2;fc), makes the right-hand 
side of (8) discontinuous and asks for a redefinition of the requirements (7). 
In this paper, we look for a locally Lipschitz radially unbounded non-negative 
functions Pk which satisfy 

Pk{zk) — and e dPk{zk) if and only if z^ & Ak, (18) 

where dPk{zk) is the Clarke generalized gradient (see Appendix B for a def- 
inition) which is needed since Pk{zk) is now not continuously differentiable. 
Similarly to (7), we are asking Ak to be the set of all local and global minima 
for Pk{zk). 

A candidate function Pk{zk) with the properties (18) and such that a function 
Xk exists for which (8), (17) hold, is the function 

P rzkj 

Pk{zk) = V / c^{s)ds, (19) 
i=i 

^The use of dynamic quantizers can lead to asymptotic results. See [10, 29] for a 
few results for discrete-time systems. Dealing with continuous-time systems and dynamic 
quantizers poses a few extra challenges, which are not addressed in this paper. 
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Figure 1: The graph of Pk{zk) with G M and A = 1. 



where Zkj is the jth component of the vector Zk G W (see Fig. 1 for a picture 

of Pk{Zk))- 

Such a function is defined on all is radially unbounded and locally 
Lipschitz. By Rademacher's theorem ([17], Chapter 3) it is differentiable 
almost everywhere. In all the points of where it is differentiable VPfe(2fc) = 
q(2fc) i.e. (8), (17) holds with Xk = M (Id : ^ is the identity function). 
Bearing in mind the definitions (4) and (16), to satisfy (18) it is necessary 
and sufficient to set a = ^. In what follows we examine the evolution of the 
system (6) under the control law: 

M 

Ui = -Y, dikqizk) , t = l,2,...,N. (20) 

k=l 

3.2 Closed-loop system 

Similarly to (10), we write the quantized control law in compact form as: 

u = -iD0lp)qiz) , (21) 
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where q{z) — {q{zi)'^ . . . ^{zm)^)^ ■ The closed-loop system then takes the 
following expression: 

i = fiO + giOi-iD^ipHz)), 

where z = {D^ (g) Ip)x and the maps /, g, h are as in (6). 



3.2.1 Control scenario and implementation 

Before proceeding to the analysis of the system, it is important to motivate 
in more detail the control scenario we consider and how the overall control 
scheme is implemented. 

For each pair of neighboring agents, one of the two is equipped with a sensor 
which continuously take the relative measurement with respect to its neigh- 
bor, e.g., a sonar or a radar. Not all the agents are equipped with these 
sensors since they might have very dedicated tasks in the formation and 
space must be saved for other hardware needed to accomphsh these tasks. 
On the other hand, since these agents need information to maintain their 
positions in the formation, they receive such information in quantized form 
from their neighbors via a digital communication channel. 

The implementation of the control law in (22) can be given by the quantization- 
based distributed control protocol as follows. 

{Initialization.) At time to = 0, all sensors measures ^fe(^o), k = 1,2, ... ,p. 
The processing units collocated with the sensors computes q{zk{to)) and the 
resulting value is broadcasted to its neighbors. Each agent i computes the 
local control law Ui as in (20) and the control value is held until new in- 
formation is available. Note that the closed-loop system evolves according 
to 

^ — ^(0 + liv<8)f 
for all t > to until new information is available. 

{Quantization-based transmission and control update.) Let i = 1 and let t^ 
be the smallest time at which t£ > and the processing unit of a sensor in 
the k-th edge detects that q{zk{t£)) ^ q{zk{ti^i)) . In this case, the quantized 
information q{zk{te)) is transmitted to its neighbor and the local control law 
of the i-th and j-th agents, where {i,j) is the pair of nodes linked by the 
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k-th edge, is updated by 

Ui{ti) = Ui{te-i) - dik{c[{zk{ti)) - q{zk{ti-i))) 

There is no other information exchange and hence the rest of the agents 
maintains their local control values. The local control law is now fixed until 
new information is transmitted again. For all t > and until this new 
transmission occurs, the evolution of the closed-loop system is given by 

The quantization-based event-triggered control update process is iterated 
with the index value £ incremented by one. 

A few remarks are in order: 

(i) The construction outhned above results in a sequence of unevenly-spaced 

sampling times ti, £ G N, at which sensors located at the systems broad- 
cast quantized information to neighboring systems. This information 
is used by local controllers to update the control value. The control 
laws turn out to be piece-wise constant functions of time whose value 
is updated whenever new information is received. 

(ii) Notice that even the agent which measures zj. implements a control law 

in which q{zk) is used instead of z^ itself. This is mainly motivated by 
our need to preserve a "symmetric" structure in the closed-loop system. 
In fact, given agents i,j and their relative distance Zk, in the case of 
unquantized information, agent i would use Zk in the control law, and 
agent j, —z^- Similarly, in the case of quantized measurements, it is 
very helpful in the analysis to employ q{zk) in the control law for agent 
i and —q{zk) in the one for agent j. 

(iii) In the quantization-based transmission and control protocol described 
above, the solution of the closed-loop system is not prevented to evolve 
along a discontinuity surface. In practice, due to delays in the transmis- 
sion and in the implementation of the control law, this could result in 
chattering, which is of course undesirable in the present context, since 
it would require fast information transmission. Nevertheless, in [13] a 
new class of hybrid quantizers have been considered which prevent the 
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occurrence of chattering. This class could also be used for the problem 
at hand in this paper, but this is not pursued further for the sake of 
brevity. 



3.2.2 A notion of solution 

The system (22) has a discontimious right-hand side due to the presence of the 
quantization functions and its analysis requires a suitable notion of solution. 
In this paper we adopt Krasowskii solutions. In fact, it was shown in [13] that 
Caratheodory solutions may not exist for agreement problems. Moreover, 
Krasowskii solutions include Caratheodory solutions and the results we derive 
for the former also holds for the latter in case they exist. 
Denoted by X{t) = F{t,X) the system (22), a function X(-) defined on an 
interval / C R is a Krasowskii solution to the system on I if it is absolutely 
continuous and satisfies the differential inclusion ([25]) 



for almost every (a.e.) t E I. The operator co(S') denotes the convex closure 
of S, i.e. the smallest closed set containing the convex hull of S. Since 
the right-hand side of (22) is locally bounded, local existence of Krasowskii 
solutions is guaranteed ([25]). 

The differential inclusion corresponding to the system (22) can be written 
explicitly. More precisely, for every k G {1, 2, . . . , M} and i € {1, 2, . . . ,p}, 
we observe that /Cq(r) is given by 



where K.q{z) :— X ^=1 /Cq(-Zfe), iCq{zk) :— X^=i /Cq(-Zfej)- Note that we have 
used the calculus rule for the set- valued map /C 



X{t) e /C(F(t, X)):=f]co {F{t, B{X, S))) 



(23) 




Using /Cq, the differential inclusion (23) for (22) can be written as 



(24) 



))] = giOi-iD ^ I,))Cci{z)) 
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(see also [19, 34, Theorem 1]). The Krasowskii solutions to (22) are also 
Filippov solutions as it follows from [25, Lemma 2.8] for a piecewise con- 
tinuous vector field F. Since every Caratheodory solutions to (22) is also a 
Krasowskii solution to (22), the stability properties of the Krasowskii solu- 
tions are also inherited by the classical Caratheodory solutions [25] in case 
the latter exist. 

3.2.3 Analysis 

RecaUing that 

and bearing in mind (11), the system (22) in the coordinates {z,$) writes as 

i = f{i)-rgm-{D®lM^)). 

Even the system above is discontinuous and again its solutions must be in- 
tended in the Krasowskii sense. It is straightforward to verify that, given any 
Krasowskii solution {x,C,) to (22), the function (z,^) = {{D^ ® Ip)x,^) is a 
Krasowskii solution to (25). The differential inclusion corresponding to (25) 
is easily understood from (24). In what follows we investigate the asymptotic 
properties of the Krasowskii solutions to (25) and infer stability properties 
of (22). 

A few notions of nonsmooth control theory which are used in the proofs are 
recalled in the Appendix B. The first fact we notice is the following: 

Lemma 1 Let Assumption 1 hold and let the communication graph G he 
undirected and connected. Then any Krasowskii solution to (25) converges 
to the set of Krasowskii equilibria: 

{{z,i) : e = 0, 0e(L>®7p)/Cq(^)}. (26) 

Proof: To analyze the system (25) we consider the Lyapunov function 

N M 
i=l k=l 

^ V ' ^ V ' 

5(0 P{z) 

k=l 3=1 
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The function is a locally Lipschitz and regular function. In fact, each term 
lo''^ <i{s)ds is convex and as such it is regular ([16, Proposition 2. 3. 6], [19]). 
Then the sums Pk{zk) and ^^^iPkizk) are also regular. The function 
V(yZ,^) is nonnegative and vanishes on the set of points such that C = 
and Zkj e [—a, a] for all A; e {1, 2, ... , M} and all j G {1, 2, . . . 
In order to apply the LaSalle's invariance principle for the differential in- 
clusions as given in Theorem 4 in Appendix B, we analyze the set- valued 
derivative V with respect to (25) as follows. Define 

V{z,i) = {a e M : 3w e ICF{z,^) s.t. a = {p,w), Vp G dV{z,^) } . 

where (■, ■) denotes the standard inner product and F{z,^) the right-hand 
side of (25). We first observe that by definition of V{z,^), dV{z,^) and 
dP{z), p G dV{z,^) implies the existence of G dP{z) such that 




Moreover, if w G )CF{z,^) then there exists Wz G K.q{z) ([25], [34]) such that 

Let now p G dV{z, ^) and w G JCF{z, ^) and write 

{p,w) = {VS{0,f{0 + 9{0i-D®Ip)^z)+ 

{p„{D^^Ip)h{0) 
A (27) 

where the inequality is a consequence of (2). Suppose now that for some 

(z,^), V{z,^) ^ 0. Then, for every a G V{z^^^ and for every p G dViz,^^ 
there exists w G }CF{z,^) such that a = {p,w). By definition of q{z) and 
P{z), dP{z) = /Cq(^) ([25], [34]). Then a = {p,w) holds in particular when 

_ f Pz \ _ f '^z \ 
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with Wz e K,<\{z). Then (27) becomes: 

N 
i=l 

N 

1=1 

Hence, for all {z,^) such that V{z,^) ^ 0, we have that 

N 

V{z, = {aeM:a<-^ Wi{^i)} . (28) 

Since |\/(^(t), ^(t)) G V{i{t),z{t)) C (-cx),0] for almost every t, V{z{t),i{t)) 
cannot increase, and any Krasowskii solution [z{t),^{t)) is bounded. Hence, 
{z{t),^{t)) exists for all t. 

Given any initial condition (2;(0), ^(0)), the set S such that V{z, ^) < V{z{{i), 
^(0)) is a strongly invariant set for (25) which contains the initial condi- 
tion. An application of Theorem 4 in Appendix B shows that any Kra- 
sowskii solution converges to the largest weakly invariant set contained in 

5n{(^,0 : oeF(z,e)}. 

Moreover, in view of (28), the set Z of points (2;,^) such that e V{z,^) 
is contained in the set of points such that ^ = 0. Hence, any point of the 
largest weakly invariant set contained in Sr[Z is such that ^ = 0. Pick a point 
{z, 0) on this invariant set. Then in order for a Krasowskii solution to (25) 
starting from this point to remain in the invariant set, it must be true that 
G /(0)+^(0)(-(L'(8)/p)/Cq(^)) = g{Q){-{D®Ip)lC(\{z)). Since the matrix 
(7(0) is full-column-rank (recall that each QiiQ) is full-column rank), the in- 
clusion above requires the existence of G /Cq(2;) such that {D®I^Wz — 0. 
In other words, the largest weakly invariant set included in iSflZ is contained 
in the set (26). Finally, observe that, taken any point in the set (26) as initial 
condition for (25), at least a Krasowskii solution (yz{t),^{t)) originating from 
this point must coincide with the trivial solution, i.e. {z{t),^{t)) = (0, 0) for 
all t. Hence, any point in (26) is a Krasowskii equilibrium for (25). □ 
It is now possible to prove the following: 
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Theorem 2 Let Assumption 1 hold and let the communication graph G he 

undirected and connected. Let v : M.^ MP be a bounded and piecewise con- 
tinuous function and A be a positive number. Then any Krasowskii solution 
to (22) converges to the set 

:e = 0, z&{Ai X ...xAm),z={D^ ®Ip)x}, (29) 

where the sets Ak 's are defined in (16), with a — A/2. Moreover, limt_^+oo[i(i) — 
lN®v{t)]^Q. 

Proof: Consider any Krasowskii solution {x{t),^{t)) to (22), whose 
existence is guaranteed locally. It can also be extended for all t G [0, +oo). 
In fact suppose by contradiction this is not true, i.e. (a;(t),^(t)) is defined on 
the interval [0,t/), with tf < +oo. Define {z{t),i{t)) = ({D'^ (g> Ip)x{t),i{t)). 
{z{t),^{t)) is a Krasowskii solution to (25). As proven before, such a solution 
is bounded on its its domain of definition. Since by (22) x{t) = h{$,(t)) +v{t) 
and both the terms on the right-hand side are bounded, then x{t) grows 
linearly in t and therefore it must be bounded on the maximal interval of 
definition, i.e. tf = +oo. Hence both {x{t),^{t)) and {z{t),^{t)) = {{D^ ® 
Ip)x{t),^{t)) are defined for all t. Moreover, by Lemma 1, z{t) = {D'^i^Ip)x{t) 
converges to the set of points (26), i.e. to 

{(x, : e = 0, e (D 8) /p)/Cq(2), z = {D^® Ip)x} . (30) 

Let (x, 0) belong to the set (30). Then z = {D^ Ip)x, i.e. z belongs to 
the span of (g) Ip and there exists Wz G A^q(-2) such that {D ® Ip)wz ~ 0. 
The two conditions imply that {wz, z) — 0. We claim that then necessarily 
2; e ^1 X ... X Ami with the sets AiS given in (16). In fact, if this is not 
true, then there must exist a pair of indices j, such that \zkj\ > a. This 
implies that the entry k + j oi the vector Wz is different from zero and also 
Wz,k+j ■ Zkj > 0. Moreover, since Wz G }Cq{z), for any pair of indices i,i such 
that i ox I ^ j, Wz,i+e ■ Zi£>Q. This contradicts that {wz, z) —Q. Then 
we have proven that the set (30) is included in the set 

{(a;, : = 0, ^ G A x . . . x Am, z ^ {D'^ ® Ip)x} . (31) 

Hence, any Krasowskii solution {x{t),^{t)) to (22) converges to a subset of 
(31). 

As for the second part of the statement, any Krasowskii solution to (22) is 
such that x{t) — lN<S>v{t) — h{$,(t)), and since we have proven that $,(t) — )> 
as t ^ 00, we have also proven that limt^+oo[2;(^) — Iat (8) v{t)] — 0. □ 
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3.3 Examples 

We provide three examples of application of the quantized agreement result 
described above. 

Agreement of single integrators by quantized measurements. We specialize 
the proof of Theorem 2 to the agreement problem for single integrators. This 
problem stems from the case when in (22) the mapping from u = — (D ® 
Ip)q{z) to y = h{^) is an identity operator. The closed-loop system (22) 
reduces to: 

X = -{D^Ip)q{z) , . 

z = iD^®Ip)x ^"^^^ 

which using the variables z becomes 

i = -{D^ (g) Ip){D (8) Ip)q{z) = -{D^D (g) Ip)q{z) (33) 

Wc analyze this system using the function P{z) introduced above. Since 
dP{z) = )Cq{z), then for all z such that P{z) ^ 0, we have 

P{z) ^ {aeR:3w, e ICqiz) s.t. a = - \\{D Ip)w^\\^}. 

Hence, arguments as in Lemma 1 give that all the Krasowskii solutions to 
(33) converge to the set of points {z : e {D <S) Ip)ICq{z)}. On the other 
hand, by Theorem 2, any Krasowskii solution x{t) to (32) is such that z{t) — 
{D'^ (g) Ip)x{t) converges to {x : G (D Ip)]Cq{z), z = {D'^ (g) Ip)x} which 
is included in the set {2; : z & Ai x . . . x Ap, z = {D^ (g) Ip)x}. 
Let X be any Krasowskii solution to (32) with z = {D^ (g) Ip)x. Take any 
two variables Xi,Xj whose agents are connected by the edge k. Consider for 
the sake of simplicity that each quantizer has the same parameter A. Then 
Zk — — Xj converges asymptotically to a square of the origin whose edge 
is not longer than A. If the agents are not connected by an edge but by 
a path, then each entry of Xi — Xj is in magnitude bounded by A • with 
d the diameter of the graph. The result can be compared with Theorem 4 
in [20]. One difference is that, while trees are considered in [20], connected 
graphs are considered here. Moreover, in [20] the scalar states are guaranteed 

to converge to a ball of radius ^j^ A. Hence, denoted by p the ratio 

A ^^-^(g^ri) considered the bound M < N — 1, any two states Xi, Xj may 
differ for 2pA^/N — 1. The passivity approach considered here yields that 
they differ for not more than d ■ A, where d grows as 0(plog(A'")) ([15]) 
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for not complete and regular graphs (graphs with all the nodes having the 
same degree), thus leading to a smaller region of convergence, the quantizer 
resolution A being the same. 

Agreement of double integrators by quantized measurements Consider the 
case of N agents modeled as 

Xi^fi, 1^1,2,. ..,N, (34) 

with Xi, fi e R^, for which we want to solve the agreement problem with 
quantized measurements. This means that all the agents should practically 
converge towards the same position and also asymptotically evolve with the 
same velocity v. The preliminary feedback ([!]) 

fi = -Ki{xi -v)+v + Ui, Ki^Kj, (35) 

with Ui to design, and the change of variables = Xi — v , makes the closed- 
loop system 

Xi = ii + v 

ii = -Ki$,i + Ui 

Vi = 

passive with storage function Si{^i) = \CT^i ^^^^ = ^^iil^i- The 

system above is in the form (1). Theorem 2 guarantees that the Krasowskii 
solutions of (34), (35), (21) converges asymptotically to the set (29) and that 
all the agents' velocities converge iov. In other words, the formation achieves 
practical position agreement and convergence to the prescribed velocity. 

Remark 2 (Consensus for double integrators with velocity feedback) 

A different but related consensus problem consists of designing local con- 
trollers in such a way that each double integrator converges asymptotically 
to the same position and velocity. In this case, no external reference ve- 
locity is provided and the velocity to which all the systems converge is the 
average of the initial velocities ([41]). The controller which guarantees this 
coordination task uses both position and velocity feedback (observe that the 
communication graphs for the position measurements and for the velocity 
measurements can be different). It then makes sense to consider the problem 
in the presence of quantized relative position and velocity measurements. 
This has been investigated in [30]. 
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The case of unknown reference velocity. If the reference velocity v is 
not available to all the agents, then [4, 6] suggest to replace it with an 
estimate which is generated by each agent on the basis of the current available 
measurements. Here we examine this control scheme when the measurements 
are quantized. We consider the special case in which the unknown reference 
velocity is constant. Then each agent i, with the exception of one which 
acts as a leader and can access the prescribed reference velocity f , use an 
estimated version of f , namely Vi has to be generated on-line starting from 
the available local measurements. The agent's dynamics (1) becomes 

Xi = yi + Vi 

ii = f^{i^) + gi{ii)Ui (36) 

Vi = hi{ii), i = 1,2, . . . , AT, 

with Vi = V Hi — 1 (without loss of generality agent 1 is taken as the leader) , 
and otherwise generated by 

Vi = KiUi 

with Aj = > and Ui as in (9). Observe that in this case, the estimated 
velocity is updated via quantized measurements. Consider the closed-loop 
system 

X — h{^)-\-lN®v 

i = fi0-9i0iD®IpHz) (37) 

V = -A{D®Ip)q{z) 

where A = diag(Ai, . . . , Ajv) and z — {D^ Ip)x. Let 

Vi{t) =v + Vi{t) - V + {vi{t) -v) =:v + Vi{t) 
where tii = 0. Rewrite the system using the coordinates z and 9 and obtain 

z = {D^®lp)[h{0+v] 

i = fiO-giOiD^ipH^) (38) 

V = A(D (g) /p)q(^) 

where in the second equation it was exploited again the fact that {D^ (8) 

Ip)lN = 0. 

One can now proceed as in the proof of Lemma 1. Consider the Lyapunov 
function 

V{z,^,v) = S{0 + P{z) + lv^A-'v 
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and let F{z, v) be the right-hand side of (38). For any p e dV{z, ^, v) and 
w e )CF{z,^,v) consider 

{p,w) = {VS{O,f{O+9m-D0lp)w,) + 

{p^,(D^^Q[hiO+v])+ (39) 
{A-^v, -A{D ® Ip)w,) 



where 



P 




and 

/ {D^®I,)[h{0+v]\ ( 
As in Lemma 1 one proves that (j),w) < — X^^^ M^i(CO therefore that 

AT 

V{z, C,v) ^ {a eR : a < -J^Wim ■ (40) 

i=l 

Hence, any Krasowskii solution {z{t),$^{t),v(t)) is bounded and exists for all 
t. Let S be the level set such that V{z,^,v) < V{z{0),^{0),v{0)) and Z the 
set of points {z,^,v) such that G 1^(2;, {;). Then any solution (z.^.v) 
converges to the largest weakly invariant subset contained in S HZ. Observe 
that Z C {{z,^,v) : ^ = 0}. Moreover, for a set in <S fl Z to be weakly 
invariant, it must be true that e ]CF{z, ^, v) with F{z, ^, v) the right-hand 
side of (38). These two facts together imply that there must exist G /Cq(^) 
such that {D®Lp)Wz = and additionally {D'^ ®Lp)v = 0. The latter implies 
that V = {In® Ip)c for some c e M. Since vi = 0, then on the largest weakly 
invariant set contained in 5" fl Z it is also true that v — 0. Hence it follows 
that any Krasowskii solution to (38) converges to the set 

{(z, e, i)) : e - 0, e (D ® Ip)ICq{z), v = 0} . (41) 

Note that each point in the set is a Krasowskii equihbria of (38). 
One can then focus on the system (37) and follow the same arguments of The- 
orem 2 to conclude that the solutions of the closed-loop system converge to 
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the set where all the systems evolve with the same velocity, achieve practical 
consensus on the position variable and the estimated velocities Vi converge 
to the true reference velocity v. 

Proposition 1 Let Assumption 1 hold and let the communication graph G 
he undirected and connected. Let v eW he a constant vector and A a positive 
numher. Then any Krasowskii solution to (37) converges to the set 

{{x,^,v) : ^ ^0,z e {Ai X . . . xAm):Z^ {D'^ ® Ip)x,v ^ 1n®v}, (42) 

where the sets Ak's are defined in (16), with a = A/2. In particular, 
\\rsit^+^[x{t) - Iat (g) w] = 0. 

Remark 3 (Velocity error feedback) Instead of the control law (10) 
u = —{D® Ip)ip{z), the control law proposed in [6] considers an additional 
velocity error injection (namely, ^^^^(j) {xj — Xi), with Af{i) the set of neigh- 
bors with respect to which the agent i can measure the relative velocity). 
This modified control law guarantees velocity tracking (with time-varying 
reference velocity) and agreement of the variables x without relying on the 
convergence of the estimated velocity to the actual value. However, the use of 
this additional velocity feedback term in the presence of quantization poses 
a few additional challenges which are not tackled in this paper. See also 
Remark 2 for more comments in this respect. 



4 Quantized synchronization of passive sys- 
tems 

We turn now our attention to the systems in (12) where the control law that 
we consider is a static quantized output-feedback control law of the form 

u = -{D (8) Ip)q(z) with z = (D'^ (g) Qw. (43) 

The overall closed-loop system is 

^ = {I^ ® A)^ - {In ® B){D ® Ip)q{z) . . 

z = {D^ (^Ip)w = {D^ (^Ip){lN(^C)^. ^ ^ 

Applications where synchronization problems under communication constraints 
and passivity are relevant are reviewed in [22]. Later in this section, we 
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briefly discuss another example where the use of quantized measurements for 
synchronization can be useful. 

To study the robustness of the synchronization algorithm to quantized 
measurements we need a more explicit characterization of the exponential 
stability of (15). To this purpose we introduce a different Lyapunov function 
which is characterized in the following lemma. As we consider time-invariant 
graphs, observability can be replaced by a detectability assumption. 



Lemma 2 Let (C, A) be detectable and U — In — 



N ■ 



The integral 



r+oo 

R:^ (n(g)7„)^e^%^"(n®7„)ds, 
Jo 



with A as in (15), is finite and satisfies 
( exp(A-A2SC)s .. 



\R\\ < 



r+oo 

Jo 



Or 



\ 



Onxn 

Moreover, the Lyapunov function 
satisfies the following: 

cm'<u{o<c,\\i\\' 

vt/(0-ie<-||eiP 

for each ^ e 



exp(A - XnBC)s J 



ds. (45) 



(46) 



Proof: The proof is given in the Appendix C. □ 
The first fact we prove about (44) is that the control law (43) achieves 
practical synchronization of the outputs: 

Proposition 2 Let Assumption 2 hold and let the communication graph G 
be undirected and connected. Then any Krasowskii solution to (44) converges 
to the largest weakly invariant subset contained in 



{^ei?'^^ : M < -,Vfc = l,2,...,M, i = l,2,...p}, 



(47) 



with z = (D^ (g) /p)(/jv <8) C)^. 
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Proof: Any Krasowskii solution to (44) satisfies the differential inclusion 

e e (/iv ® A)i - {In ® ® Ip)ICq{z). 

Consider the Lyapunov function V{^) = ^^(/at^P)^. Then, for any ^ G i?^" 
and any 1/ e ]Cq{z), with z — {D^ (8) /p)(/jv <S) C)^, we have 

V^(0 VF(0-[(/iv<8^)e-(/iv<8S)(i?®/p)i^] 
= 2e^(7Ar (8) PA)C - 2^^ {In ® PS) {D (g) 7p)i/. 

Using Assumption 2 and the definition of z we further obtain that for all 
u e /Cq(2;), 

V{0 < -2e{lN®C^){D®I,)y 

= -2^^i/<0. ^ ^ 

This shows that V{^{t)) cannot increase and that ^{t) is bounded. Moreover, 
by LaSalle's invariance principle for differential inclusion (Appendix B, The- 
orem 4), any Krasowskii solution converges to the largest weakly invariant 
subset contained in 

e R^" : 3 1/ e /Cq(z) s.t. VV{i) ■ [{In ® A)i - {In ® B){D ® Ip)v\ = 0}. 

In view of (48) , any point in this set is such that Zkj Vkj = for all k = 
1,2,..., AI and for all j = 1,2, ... ,p. Since u^j G )Cq{zkj), then Zkj I'kj = 
imphes that \zkj\ < ^. This ends the proof. □ 

Remark 4 (Practical output synchronization) Similarly to the consen- 
sus problem under quantized measurements (see Section 3.3), a consequence 
of the previous statement is that any two outputs Wi,Wj practically asymp- 
totically synchronize. Namely, considered any Krasowskii solution ^{t) and 
the corresponding output w{t) = {In C)^(t), for each i = 1,2, . . . ,n and 
each t > 0, the difference \wi£{t) — Wjf{t)\ is upper bounded by a quantity 
which asymptotically converges to d^, with d the diameter of the graph. 

The proof of the proposition above clearly does not rely on the lin- 
earity of the systems but rather on the passivity property. Hence, if one 
considers nonlinear passive systems, that is systems for which a positive 
definite continuously differentiable storage function Vi{^i) exists such that 
VV^ • fi{Ci,Ui) < wj • Ui, with Wi = hi{^i), then for the closed-loop system 

— fi{^ii'^i)-i with Ui given in (43), and i — 1,2, . . . , N , it is still true 
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that the overall storage function V{^) = X]i=i ^(^i) satisfies the inequality 
V{^) < -2z'^v for all z ^ {D'^ ® IpMO and all u e JCq{z). Hence, the 
following holds: 

Proposition 3 Let Assumption 1 hold and let the communication graph G 

he undirected and connected. Then any Krasowskii solution to the systems 
(1) in closed-loop with u = —{D (g) Ip)q{z) and z = {D^ (g) Ip)h{C,) converges 
to the largest weakly invariant subset contained in the set (47). 

The next lemma states a property of the average of the solutions to (44) 
which helps to better characterize the region where the solutions converge. 

Lemma 3 Let Assumption 2 hold and let the communication graph G be 
undirected and connected. Any Krasowskii solution ^{t) to (44) satisfies 



Bearing in mind that for matrices F e ]^"»x" Q g M.p^i^ the following 
property of the Kronecker product holds: 



where in the equality before the last one it was exploited the fact that 
iJfD — O]^, which holds by definition of the incidence matrix D. Hence, 
any Krasowskii solution ^{t) is such that the average (1^ In)^{t) satisfies 




one can further show that 



dt 



e {ll^Q{lN^A)at)-B{l%D®Ip))Cq{z) 



(49) 



(l?;®4){(t)=e^*(l^®/„)C(0). 



□ 
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The following result provides an estimate of the region where the solu- 
tions converge and shows practical synchronization under quantized relative 
measurements: 

Theorem 3 Let Assumption 2 hold and let the communication graph G be 
undirected and connected. Assume that (C, A) is detectable. Then for any 
Krasowskii solution ^{t) to 

i = {In ^ A)^ - {In ® B){D ® Ip)q{{D'^ ^ Ip){lN ® C)0 (50) 
there exists a finite time T such that ^{t) satisfies 



VpM 



m - (1 



N 



N 



<2^J\R\\\\B\\\\D ^ Ip\\A (51) 



for allt > T , where Ci, C2, ||-R| | are defined m (45), (46). Moreover, ^^§^^{t) 
^o{t) where C,o{t) is the solution of ^o{t) = A^Q{t) with the initial condition 
eo(0) = ^e(0). 

Proof: By definition, any Krasowskii solution ^ to (50) is such that 



^ = (n (g) /„){, with U = In - satisfies 

i e {In (8) A)i- {In ^ B){D ^ Ip)^^.{{D'^ ® ® C)0, 

where similar manipulations as in (49) were used. Moreover, any u e /Cq((i?-^§ 

Ip){lN ® C)^) is such that ||z/ - (D^ (g) Ip){lN (S) C)^\\ < \/pM^. Under 

the assumption on the detectability of (C, A), we can consider the Lya- 
punov function U{^) introduced in Lemma 2. For any ^ and any v e 
/Cq((L'^(8)7p)(7Ar®C)0, 

Vf/(0 [{In ® A)l - {In ® B){D ® Ip)v] 
= VU{1) [{In ®A)- {In B) {DD'^ ® Ip) {In ® C)]l+ 
VU{i){lN ^B){D^ Ip)[{D'^ (8) Ip){lN (8) C)i - u] 

< -[[i[[{[[i[[-[[R[[[[B[[[[D®Ip[[^A). 



T 



Hence, for ||||| > [[B[[ [[D ^ Ip[[y/j^A, 

Vt/(0[(/iv ® A)^ - {In (gB){D(g Ip)iy] < < -^U{0. 
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It follows that any Krasowskii solution converges in finite time to the set of 
points ^ such that 



iieii< 



2y^||i?|| \\B\\ ||D®/p||VpMA 



from which the thesis is proven by definition of ^. 

The proof of the final claim follows from the fact that by Lemma 3, for 
all t > 0, 

□ 

Remark 5 (Role of ||-R||) In the case A = 0, B — C — 1, the bound on R 

reduces to ||i?|| < where A2 is the algebraic connectivity of the graph. In 
this case, the size of the region of convergence in (51) resembles the estimate 
given in Theorem 1 and Corollary 1 in [13] for quantized consensus of single 
integrators. Theorem 3 can be viewed as the extension of the results in [13] 
to the problem of synchronization of linear multi- variable passive systems by 
quantized output feedback. 



4.1 Examples 

In the following examples, we discuss how synchronization with quantized 
measurements can play a role in a decentralized output regulation problem 
in which heterogeneous systems asymptotically agree on the trajectory to 
track. 

Output synchronization for heterogeneous linear systems. In [43] (see also 
[3] , Section 3.6) the following problem is investigated. Given N heterogeneous 
linear systems 

— FiXi -\- GiUi 
Vi = HiXi, i = l,2,...,N ^ ' 

with (Fj, Gi) stabilizable and {Hi, Fi) detectable, and a graph G (which here, 
as usual in this paper, we assume static undirected and connected), find a 
feedback control law Ui for each system i (i) which uses relative measurements 
concerning only the systems which are connected to the system i via the graph 
G and (ii) such that output synchronization is achieved, i.e. limt^oo WUiif) — 
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y,(i)|| = OforalH,ie{l,2,...,7V}. 

Excluding the trivial case in which the closed-loop system has an attractive 
set of equilibria where the outputs are all zero, the authors of [43] show that 
the output synchronization problem for N heterogeneous systems is solvable 
if and only if there exist matrices S, R such that limt_^oo I IViit) ~ R^'^^WqW — 
for each i e {1, 2, . . . , N}, for some Wq. Moreover, provided that a{S) C jR, 
the controllers which solve the regulation problem are 

Xi — FiXi + GiUi + Li{yi — CiXi) 

yi = HiXi (53) 

Ui = Ki{xi - Ui^i) + Fi^i 

where e MP are the exosystem states that synchronize via communication 
channels and are described by 



i = {In ^ S)^ - {In ^ B){D ® Ip)z 
z = {D'^ ®Ip){lN®C)^, 



T^j.n ^n,c (54) 



where D is the incidence matrix associated to the graph, the pair (C, S) is 
detectable the matrices Lj, Ki are such that Fi + GiKi, Fi + LiHi are Hurwitz, 
and rij, Fj are matrices which solve the regulator equations 

F.Iii + GiVi = HiS 
HiUi = R. 

The controllers (53)-(54) are a modified form of the ones in [43, Eq. (10)] 
where in the latter, the local controller communicates the entire exosystem 
state to its connecting nodes. When the relative measurement Zk is trans- 
mitted via a digital communication line, then this information is quantized 
and the variable z in the controller (53)-(54) is replaced by its quantized 
form q(^). 

Let the eigenvalues of S have in addition multiplicity of one in the minimal 
polynomial, so that we can restrict S to be skew-symmetric without loss of 
generahty and B — C^. Then the exosy stems 

Wi = C^i z = l,2,...,Ar ^ > 

trivially satisfy Assumption 2. Then Theorem 3 applies and the sohitions ^j, 
i = 1, 2, . . . , A^, of (54) practically synchronize under the quantization of z. It 
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is then possible to see that the closed-loop system of (52) and the controllers 
(53)-(54) with z replaced by o^z) achieves practical output synchronization. 
This follows from similar arguments as in [43, Theorem 5] where [43, Theorem 
1], which is used in the proof the theorem, is replaced by Theorem 3. 

Before ending the section, we remark that Theorem 3 also holds under a 
slightly different set of conditions which do not require passivity. 

Assumption 3 Let (A, 5, C) he stabilizable and detectable, and assume that 

[7p + A^G][7^ + A2G]-i (56) 

is strictly positive real where G{s) = C{sl — Ay^B is the transfer function 
of (12) and Xn is the largest eigenvalue of L. 

Under Assumption 3, the results in Theorem 3 still hold mutatis mutandis. 
Indeed, by the multivariable circle criterion in [27, Theorem 3.4], (A — XiBC) 
is Hurwitz for every non-zero eigenvalue Aj of L. This implies that (15) is 
exponentially stable (this is evident from the proof of Lemma 2 - see the 
Appendix) and Lemma 2 and 3 continue to hold. As a consequence the 
proof of Theorem 3 holds word by word under the assumption that {A, B, C) 
is minimal and Assumption 3 holds. 

The case of output synchronization with filtered and quantized signals. 
As a concrete example to the case of exosystems satisfying Assumption 3, 
we consider again the closed-loop systems in the previous example where 
the heterogenous linear systems (52) are interconnected with the controllers 
(53)-(54) with 








u 







' " 


s = 


—OJ 








,B = 


1 







a 


—a 








,C= [ 1 ] 



(57) 



The system (S, B, C) can be considered as a cascade interconnection of a 
second-order oscillator with frequency u and a low-pass filter with a cut-off 
frequency a, and its transfer function is given by 



G{s) 



as 



(s2+a;2)(s + a) 



Using the above {S,B,C), the interconnected exosystems (54) with quan- 
tized measurement q(^) resemble a network of oscillators where the relative 
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measurements Zk are filtered and quantized. In the limiting case a — > oo, the 
exosystems are given by (55) where 



A 



oj 
-u 




1 



,C=[0 1]; 



(58) 



and it satisfies Assumption 2. A direct application of Theorem 3 shows that 
(51) holds with 



\R\\ < 



( 



exp 



V 



u 
-u -A2 

On y n 



exp 



0. 



uj 



\ 



J 



ds. 



In particular, if A2 > 4w^, then ||i?|| < / 

On the other hand, if < a < 00, i.e., when the low-pass filter is used, then 
it can be checked that 



inf Re 

V 



1 + \j~^G{iv) 



> 



1 + \2Gitiy) 

^mi{a(jj^ - av'^f + ((cj^ + \Na)v - i^^)((cu^ + \2a)u - u^) > 0. 

Note that for a sufficiently large a > 0, the above condition holds. Thus, the 
cut-off frequency a can be designed based only on the knowledge of A2,Ajv 
and u, such that the exosystems (55) satisfy Assumption 3. 
In both cases, practical output synchronization of the closed-loop systems 
(52)-(54) with quantized q{z) is obtained. 



5 Conclusions 

The passivity approach to coordinated control problems presents several in- 
teresting features such as for instance the possibility to deal with agents 
which have complex and high- dimensional dynamics. In this paper we have 
shown how it also lends itself to take into account the presence of quantized 
measurements. Using the passivity framework along with appropriate tools 
from nonsmooth control theory and differential inclusions, we have shown 
that many of the results of [1, 38] continue to hold in an appropriate sense in 
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the presence of quantized information. We believe that the results presented 
in the paper are a promising addition to the existing literature on continuous- 
time consensus and coordinated control under quantization ([20, 13, 30, 23]). 
Many additional aspects deserve attention in future work on the topic. The 
approach to quantized coordinated control pursued in this paper appears to 
be suitable to tackle more complex formation control problems such as those 
considered e.g. in Section II. C of [1], [20], Section 4 and [41]. These possible 
extensions can also benefit from the results of [6]. 

In the paper it was not discussed whether or not the use of quantized measure- 
ments yields sliding modes. Sliding modes were shown to occur in problems 
of quantized consensus for single integrators ([13]) and hysteretic quantizers 
were introduced to overcome the problem. A similar device could prove use- 
ful in quantized coordination problems. 

The literature on synchronization and coordination problems which exploit 
passivity is rich (see e.g. [35, 38, 14, 42] and references therein) and the prob- 
lems presented there could be reconsidered in the presence of quantized mea- 
surements. The book [3] provides many other results of cooperative control 
within the passivity approach. These results are all potentially extendible to 
the case in which quantized measurements are in use. 

Acknowledgement The authors would like to thank Paolo Frasca for a 
remark on the first example in Section 3.3. 
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A Notation 

The Kronecker product of the matrices A e R"*^", B e R^^* is the matrix 



2010. 



/ On B 




A®B 



\ «ml B 




See e.g. [1, 38] for some basic properties 
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B Nonsmooth control theory tools 



A few tools of nonsmooth control theory which are used throughout the 
paper are recalled in this appendix (see [2, 19] for more details). Consider 
the differential inclusion 

X e F(x), (59) 

with F : M" — )■ 2^^" a set-valued map. We assume for F the standard 
assumptions for which existence of solutions is guaranteed ([25]). Xq G is 
a Krasowskii equilibrium for (59) if the function x(t) — Xq is db Krasowskii 
solution to (59) starting from the initial condition Xq, namely if G F{xq). 
A set S is weakly (strongly) invariant for (59) if for any initial condition 
X G 5 at least one (all the) Krasowskii solution x{t) starting from x belongs 
(belong) to S for all t in the domain of definition of x{t). Let 1/ : M" — )■ M be 
a locally Lipschitz function. Then by Rademacher's theorem the gradient of 
V exists almost everywhere. Let N be the set of measure zero where W{x) 
does not exist. Then the Clarke generalized gradient of F at x is the set 
dV{x) — co{limj_^+oo '^V{xi) : — >■ x, Xi ^ S , Xi ^ N} where S is any set 
of measure zero in M". We define the set- valued derivative of y at a; with 
respect to (59) the set V{x) = {a G M : 3f G ICf{x) s.t. a = p ■ v, Vp G 
dV{x)}. The definition of regular functions used in the following nonsmooth 
LaSalle invariance principle can be found e.g. in [2]: 

Theorem 4 ([2, 18]) Let V" : M" — )■ M &e a locally Lipschitz and regular 
function. Let x E S, with S compact and strongly invariant for (59). Assume 

that for all X E S either V{x) = % or V{x) C (— oo, 0]. Then any Krasowskii 
solution to (59) starting from x converges to the largest weakly invariant 

subset contained m <S fl {x G : Q &V{x)}, with the null vector in R". 

C Proof of Lemma 2 

Proof: Following [33] , Theorem 3 (see also [21]) , we introduce the NxN 
nonsingular matrices ^ 

T={1n/VNv2...Vn),T-'={1n/^ W2 ... WNf 

^The matrices T, transform the Laplacian matrix L = DD^ into its diagonal form. 
The columns of T form an orthonormal basis of M^. 
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and notice the following: 



( 0. 



On> 



0. 



exp(yl - \2BC)s 



Onxn 



0. 



(T-i®4)(n®4) 



(60) 



exp(yl - XnBC)s J 



where < A2 < ... < Aat arc the non-zero eigenvalues of DD^. Since 
{A, B, C) is passive and (C, A) is detectable, then the matrices A — XiBC, 
i = 2, . . . , N are Hurwitz. This implies that only exponentially stable modes 
are present in e'^*(n ® /n), and therefore the integral which defines R exists 
and is finite. 

Using the transformation matrix T, a routine computation shows that 

(n®4)^e^%^^(n®4) 





^ Onxn 
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Ojjxn ^ 
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Taking the norm of the matrix, 

||(n®4)^V^e^^(n®4)|| < 



/ exp{A-X2BC)s ... Onxn \ 

\ Onxn ■ ■ ■ exp{A - XnBC)s J 



from which (45) follows. 
Rewrite the function U{^) as 

/ |^(n 4)^e^^(--*)e^("-*) (n ® Qidr 

/+0O 

where ^(r; ^, t) is the solution to (15) at time r starting from the initial condi- 
tion ^ at time t. Following standard converse Lyapunov theorem arguments 
(see e.g. Khalil, Theorem 4.12) one easily proves that 



Moreover, 



ci\\i\\'<u{i)<c,m 



vuiOA^ = [e^(n®4)^e^^^e^^(n®4)e 

= -iieip- 



s=+oo 



J s=0 



□ 
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